Testing the Non-universal Z' Model in B s — >■ <pK Decay 



JuanHua 1 f C. S. Kim 2t , YingLi 1 ' 2 * 

1. Department of Physics, Yantai University, Yantai 264-005, China 
2. Department of Physics and IPAP, Yonsei University, Seoul 120-479, Korea 

February 12, 2010 

Abstract 

The branching ratio and direct CP asymmetry of the decay mode B s — > <j>7t° have been calculated within the QCD 
factorization approach in both the standard model (SM) and the non-universal Z' model. In the standard model, the 
CP averaged branching ratio is about 1.3 x 10~ 7 . Considering the effect of Z' boson, we found the branching ratio 
can be enlarged three times or decreased to one third within the allowed parameter spaces. Furthermore, the direct CP 
asymmetry could reach 55% with a light Z' boson and suitable CKM phase, compared to 25% predicted in the SM. 
The enhancement of both branching ratio and CP asymmetry cannot be realized at the same parameter spaces, thus, 
if this decay mode is measured in the upcoming LHC-b experiment and/or Super B-factories, the peculiar deviation 
from the SM may provide a signal of the non-universal Z' model, which can be used to constrain the mass of Z' boson 
in turn. 
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Although most of the experimental data are consistent with the standard model (SM) predictions, it is believed 
that the SM is just an effective theory of a more fundamental one yet to be discovered. One way of searching for new 
physics beyond the SM is by studying the rare B decay modes, which are induced by flavor changing neutral current 
(FCNC) transitions, since such rare decays arise only from the loop level within the SM. Over the years, many studies 
have been made to predict the branching ratios and CP asymmetries of B decays in the SM and in new physics (NP) 
models, such as supersymmetry and etc. Although the presence of NP in the b sector is not yet firmly established, there 
exist several signals which will be verified in the forthcoming LHC-b experiment and super-B factories. Therefore, it 
is interesting to explore as many rare decays as possible to find an indication of NP. 

Additional U(l)' gauge symmetries and associated Z' gauge bosons \\'\ could appear in several well motivated 
extensions of the SM. Searching for an extra Z' boson is an important mission in the experimental programs of Tevatron 
and LHC. One of the simple extensions beyond the SM is the family non-universal Z' model, which could be naturally 
derived in certain string constructions J2), E6 models J3) and so on. It is interesting to note that the non-universal Z' 
couplings could lead to FCNC in the tree level as well as introduce new weak phases, which are essential in inducing 
the CP asymmetries. The effects of Z' in B sector have been investigated in a number of papers, such as Refs. (4][5). 
The recent review about Z' in detail is referred to Ref . J6l . 

In this work, we will address the effect of the Z' in the rare decay mode B s — > It is expected to have a small 
branching ratio in the SM because it is an electro-weak penguin dominated process and mediated by b — > sqq. In 
dealing with the two body charmless non-leptonic B decays, many approaches have been proposed, such as the naive 
factorization, the QCD factorization (QCDF) approach [SJ , the perturbative QCD (PQCD) approach and the soft 
collinear effective theory (SCET). In previous studies, the branching ratio is shown to be about 10~ 7 in the SM, both 
in the QCDF approach O and in the PQCD approach (9). For completeness, we would first calculate the mode within 
the SM, before discussing the effect of the new physics. Since there is no annihilation contribution in this decay, we 
will adopt the QCDF approach. 

We start from the relevant effective Hamiltonian given by: 



V ub V* us [C l O p 1 + C 2 O p 2 - V tb V* £ QOi 



(1) 



The explicit form of the operators 0\ and the corresponding Wilson coefficients Q at the scale of pi = nif, can be found 
in Ref. |[T0l . V u u\ b , V u u\ s are the Cabibbo-Kabayashi-Maskawa (CKM) matrix elements. 

In the QCDF approach, the contribution of the non-perturbative sector is dominated by the form factors of B s —> 
transition and the non-factorizable impact in the hadronic matrix elements is controlled by hard gluon exchange. The 
hadronic matrix elements of the decay can be written as 

l d£, f dx f ~ dyif(£,x,y)<bB(S)*t(x)* x (y). (2) 
Jo Jo 

Here T?, and 7] denote the perturbative short-distance interactions and can be calculated perturbatively. <!>x(x) (X = 
B s ,7t,(j>) are the universal and non-perturbative light-cone distribution amplitudes, which can be estimated by the light 
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cone QCD sum rules. Following the standard procedure of QCD factorization approach, we can write the decay 
amplitude as 

= ^ E I V p6 V>f(/i)(0»°|O i |fi>. (3) 

V ^ p=u,c i 

where {M\M2\Oi\B)p is the factorizable matrix element, which can be factorized into a form factor times a decay 
constant, and the coefficients a,- (z = 1 to 10) can be found in Refs. |7][8)- Note that in dealing with the hard-scattering 
spectator interactions in the QCDF, there is an infrared endpoint singularity, which can only be estimated in a model- 
dependent way with a large uncertainty. In Refs. QUI, this contribution is parameterized by one complex quantity 

X H =(l+p H e i *»)ln'^, (4) 
A/, 

where A/, = 0.5 GeV, 0# is a free strong phase in the range [—180°, 180°], and p# is a real parameter varying within 
[0,1]. 

Finally the decay amplitude can be given as 

V u i,V* s (a2[u] +a^[u] — a^[d] — as[u] + as[d] — aj[u] — —aj[d] + ag[u] + —ag[d]) 



■ V c i,V* s (03 [u] — «3 [d] — a 5 [u] + «5 [d] — a-j [u] — — a-/ [d] + ag [u] + — ag [d] ) 



(5) 



where the symbols u and d in square brackets indicate the component of the meson n°. In the SM, a,-[«] = a/[d] = a,- 
therefore, we get the simplified formula for the decay amplitude: 



° <j)7Z ) = -i^LmlfnA^iO) x \v ub V: s (a 2 - \a n + la 9 ) +V cb V*(-l<z 7 + |a 9 ) 



(6) 



after utilizing 2m* (e* ■ pb) ~ rrig. The branching ratio takes the form 

where Tg is the B s meson lifetimes, and \P C \ is the absolute value of two final-state hadrons' momentum in the B s rest 
frame. We can also define the direct CP asymmetry as: 

CP \£/(B0^<p7z°)\ 2 + W(B^<i>n )\ 2 ' 
Note that in the naive factorization there is no CP asymmetry because of none existence of any strong phase, which is 
a key factor in producing a direct CP asymmetry. 

For the numerical calculation, with the input parameters listed in Table. Q] the averaged branching ratio and direct 
CP asymmetry of decay B s — > <j)K° obtained in the SM are 

38(B S — >• (j>n°) = 1.3 xl0~ 7 , 
A CP (B S <j)7t°) = 25%, (9) 
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which have not yet been measured in the Tevatron experiments. However, the order of magnitudes should be measured 
easily in the LHC-b experiment and/or Super B-factories in future. Because we used the updated parameters, the 
branching ratio is slightly larger than that predicted in Ref. (8), and the CP asymmetry agrees with each other. The 
results also agree with the predictions from the PQCD J9J as well. Here we will not tend to discuss the uncertainties 
in our calculation, since this part has been presented explicitly in j8||. 

Now we turn to the effects due to an extra U(l)' gauge boson Z'. We start from the interactions with the new Z' 
gauge particle ignoring the mixing between Zq and Z'. Following the convention in Ref. JT|, we write the couplings 
of the Z' -boson to fermions as 

4 =g%Vit[e7 L PL + £7 R P R ]Vi, (10) 

i 

where i is the family index and \j/ labels the fermions and Plr = (1 qp According to certain string constructions 

ifTTl or GUT models lfl2l . it is possible to have family non-universal Z' couplings. That is, even though £■' are 
diagonal, the couplings are not family universal. After rotating to the physical basis, FCNC's generally appear at tree 
level in both left handed and right handed sectors, explicitly, as 

B^=V VL e^Vl, BV« =V w eV*Vl R . (11) 

For simplicity, we assume that the right-handed couplings are flavor-diagonal and neglect Bf b , thus the Z' part of the 
effective Hamiltonian for b — > sqq (q = u,d) transitions has the form as: 

JHg f = ^(^fB^y^iB^qqh-A+Bfj^y^) +h.C. , (12) 

where g\=e/ (sin 9w cos 9w ) and M z > is the new gauge boson mass. Compared with the operators existed in the SM, 
Eq. (TT~2b can be modified as 

Mff f = -^V, fe V*£(AC 3 Of +AC 5 0* + AC 7 0* + AC 9 0*) +h.c. , (13) 
where O q { (i = 3,5,7,9) are the effective operators in the SM, and AC, the modifications to the corresponding SM 



Table 1 : Summary of input parameters (8) 
X A p f) A^T 4) t^o X B a e a s 
0.225 0.818 0.141 0.348 250MeV 1.46ps 0.35 1/132 0.214 
fs s m Bs f n y A S ^' 

236MeV 5.36GeV 131MeV 221MeV 175MeV l.OlGeV 70° 0.34 
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Wilson coefficients caused by Z' boson, which are expressed as 



AC 3 = 

AC 5 = 

AC 7 = 

AC 9 = 



^ ( ^ )2 ^ (B "" +2BL) ' 
4 , g'M z , 2 

3y,f,y,:. [ gl M z , 

4 



' B dA B uu - B dd) I 



( 5 »/ ) 2 B ^ ~ B dd); 



(14) 

3^V M giM z / 

in terms of the model parameters at the Mw scale. While we can have Z' contributions to the QCD penguins as well as 
the EW penguins, in view of the results evaluated by Buras et. al |[T3l . we set Buf = —2Bbj, so that new physics is 
manifest in the EW penguins. Without loss of generality, we always assume that the diagonal elements of the effective 
coupling matrices Bq f are real due to the hermiticity of the effective Hamiltonian. However, there still is a new weak 
phase <p in the off-diagonal one of B^ b . The resulting Z' contributions to the Wilson coefficients are: 



AC 3j5 0, 
AC 9 7 = 4 



\v t bK\ 



--L,R-iA 



with 



v tb v* 

g'M z N 
giM z 



r>L r>L R 
B sb B dd 



v th v* 



(15) 



(16) 



To address the effect of Z' boson, we have to know the values of the AC7 and AC9 or equivalently B^ b and Bzf. 
Generally, we always expect g'/gi ~ 1, if both the U(l) gauge groups have the same origin from some grand unified 
theories. And Mz/M z i ~ 0.1 for TeV scale neutral Z' boson, which yields y ~ 10~ 2 . In the first paper of Ref. 101 
assuming a small mixing between Z Z' bosons the value of y is taken as y ~ 10~ 3 . In order to explain the mass 
difference of B s —B s mixing, we need \B^ b \ ~ IV^V^-I- Similarly, the CP asymmetry anomaly in B — > (j)K,nK can 



be resolved if 



|V t6 V,*|, which indicates 

} dd\ 



1. Above issues have been discussed widely in Ref. 0. 



Because we expect that \Bh d \ and \B^ S \ should have the same order of magnitude, we simply assume that 



III = 1111 = \ti\ = \\tf\ = e(io- 3 ,K 



-2 



(17) 



since the major objective of our work is searching for new physics signal, rather than producing acute numerical 
results. Due to renormalization group (RG) evolution from the Mw scale to m b scale, the other Wilson coefficients 
also receive the contribution of Z', however, the RG running from the m z i to Mw scale has been neglected in this work. 



The Wilson coefficients at m b and \/ Aj x m b scale have been presented in Table. 2. 



Once obtaining the values of the Wilson coefficients at the scale m b and ^/7^^m b , we can get the decay amplitude 
from the Z', analogous to Eq. ((5), as: 



A*/(B° s (pB)^H£,Pt)n°(p2)) 



-i^2m p f K (e*-p B )A B ^(0) 



V u bV* s (A«2 [u] + A«3 [u] — A«3 [d] — A«5 [u] + Aa$ [d] — Aa-j [u] — ijAa-i [d] 



- A«9 [u] 



■ —Aay [d] 



' ^ 'cbVc S {AaT,[u\ — Afl3[c/] — a$[u] + Aa^[d\ — Ag^m] — — A«7[c/] + Aagfj*] + — A«9[rf]) 



(18) 
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Figure 1 : After setting £ =0.01, the variation of the CP averaged branching ratio (left panel) and direct CP asymmetry 
(in %) (right panel) as a function of the new weak phase 0. We varied the unitary angle y G (50°, 1 10°). The horizontal 
lines are predicted in the SM. 




Figure 2: When setting y = 70°, the variation of direct CP asymmetry with the new weak phase <j), where the solid, 
dot-dashed and dashed lines correspond to £, = 0.001,0.005 and 0.01. 
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Table 2: The Wilson coefficients Q within the SM and with the contribution from Z' boson included in NDR scheme 
at the scale jx = and ji/, = VA/,m^. 



Wilson 


M = m b 




= V A h"V, 


coefficients 


C SM 


ACf 




ACf 


Ci 


1.075 


-0.006£ L 


1.166 


-0.008^ 


c 2 


-0.170 


-0.009§ L 


-0.336 


-0.014<§ L 




0.013 


0.05£ L -0.0l£ s 


0.025 


0.11^-0.02^ 


Ca 


-0.033 


-0.13^+0.01^ 


-0.057 


-0.24£ L + 0.02£ fl 


c 5 


0.008 


0.03£ L + 0.01<^ 


0.011 


0.03^ L +0.02^ 




-0.038 


-0.15£ L + 0.0l£ R 


-0.076 


-0.32^ L + 0.04^ 


Ci/a em 


-0.015 


4.18£ L -473£ R 


-0.034 


5.7^-459^ 


Q/ <Xem 


0.045 


1.18^ L -166^ 


0.089 


3.2^-355$* 


Cg/a em 


-1.119 


-561^ L +4.52^ 


-1.228 


-611§ L + 6.7<^ 


Clo/ %m 


0.190 


118£ L -0.5£ R 


0.356 


207^-1.4^ 



To study the effect of the Z' boson, by setting ^ = 0.01 and varying y within 50° to 1 10°, one can get the variation of 
the CP averaged branching ratio and the direct CP asymmetry as a function of the new weak phase <p, as shown in Fig. 
Q] where the horizontal lines are the values predicted in the SM. From these figures, we find that the branching ratio 
may become three times of that predicted in the SM or drop to one third of the SM value within the allowed parameter 
space. Moreover, as we mentioned before, we have introduced one new weak phase from the off-diagonal element of 
B^ b , which plays a major role in changing the direct CP asymmetry. The direct CP violation can reach 55% if y = 50° 
and = 70°. This remarkable enhancement will be an important signal in testing the model. Taking y= 70°, we plot 
the variation of direct CP asymmetry as a function of the new weak phase (j) with different £, = 0.001,0.005,0.01, as 
shown in Fig. [2] According to this figure, we note that the new physics effect cannot be detected if ^ < 0.001, namely 
a heavier Z' boson. If there exists a light Z' boson, the observation of this mode will in turn help us constraint the 
mass of Z'. In Fig. [3] when leaving the £, and as free parameters, and setting y = 70°, we present the correlations 
between the averaged branching ratio, direct CP asymmetry and the parameter values by the three-dimensional scatter 
plots. As illustrated in Fig. [3] the enhancement of both branching ratio and CP asymmetry cannot be fulfilled at the 
same parameter values. 

To conclude, we have calculated the branching ratio and direct CP asymmetry of the decay mode B s — >• <j)K° within 
the QCD factorization approach in both the SM and the non-universal Z' model. This approach is suitable as the decay 
mode has no pollution from annihilation diagrams. Upon calculation, we found the branching ratio may be enlarged 
three times or decreased to one third by the effect of Z' boson within the allowed parameter space. Furthermore, as the 
direct CP asymmetry is concerned, it can reach 55% with a light Z' boson and suitable CKM phase. Also, we note the 
enhancement of both branching ratio and CP asymmetry cannot be accomplished at the same parameter space. Thus, 
if this mode could be measured in the upcoming LHC-b experiment and/or Super B-factories it will provide a signal 
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Figure 3: The variation of the CP averaged branching ratio (left panel) and the direct CP violation (right panel) with ^ 
(in units of 10~ 3 ) and the new weak phase 0. 

of the non-universal Z' model, and can be used to constrain the mass of the Z' boson in turn. 
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